Microscopic Theory of Heterogeneity and Non-Exponential Relaxations in 

Supercooled Liquids 
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Recent experiments and computer simulations show that supercooled liquids around the glass 
transition temperature are "dynamically heterogeneous" Q. Such heterogeneity is expected from 
the random first order transition theory of the glass transition. Using a microscopic approach based 
on this theory, we derive a relation between the departure from Debye relaxation as characterized 
by the f3 value of a stretched exponential response function 4>{t) = ^-(^/'^kww) ^ a,nd the fragility of 
the liquid. The /3 value is also predicted to depend on temperature and to vanish as the ideal glass 
transition is approached at the Kauzmann temperature. 
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The striking universality of relaxation dynamics in su- 
percooled liquids has remained intriguing for decades. In 
addition to the overall dramatic slowing of transport as 
the glass transition is approached, one finds the emer- 
gence of a strongly non-exponential approach to equilib- 
rium when a supercooled liquid is perturbed. This con- 
trasts with the behavior of chemically simple liquids at 
higher temperatures, where only a single time scale for 
a highly exponential structural relaxation is usually en- 
countered for times beyond the vibrational time scales. 
Both the range of time scales and the median magnitude 
of the relaxation time require explanation. 

Several theoretical threads lead to the notion that 
the universal behavior of supercooled liquids arises from 
proximity to an underlying random first order transition 
[1I~H which is found in mean field theories of spin glass 
without reflection symmetry |7|-|9[|, and in mode coupling 
and density functional p^-|l^ approaches to the 
structural glass transition. This picture explains both 
the breakdown of simple collisional theories of transport 
that apply to high temperature liquids at a characteristic 
temperature Ta and the impending entropy crisis of su- 
percooled liquids first discovered by Simon and brilliantly 
emphasized by Kauzmann ]l5[ | at the temperature Tk- 
Furthermore the scenario suggests that the finite range 
of the underlying force modifies mean field behavior be- 
tween Ta and Tr- in a way that leads to an intricate "mo- 
saic" structure of a glassy fluid in which mesoscopic local 
regions are individually each in an aperiodic minimum 
but are separated by more mobile domain walls which are 
strained and quite far from local minima structures [ p^ . 
Relaxation of the elements of the mosaic, reconfiguring 
to other low energy structures, leads to the slow relax- 
ation and a scaling treatment of the median relaxation 

time yields the venerable Vogel-Fulcher law, r = roe^^^o 
[^. D, called the fragility, determines the apparent size 
of deviations from an Arrhenius law. Most recently a 
microscopic calculation of the coefficient D has yielded 



good agreement with experiments for a range of liquids 
Jirt . In this paper, we address the predictions of the 
mosaic picture for the dispersion of relaxation times. 

Relaxation in supercooled liquids is well approximated 
by the stretched exponential or Kohlrausch- Williams- 
Watts (KWW) formula = g-Ct/^^ww)", a study 
conducted by Bohmer et al. on over sixty glass for- 
mers around Tg shows /? and D are strongly correlated. 
The smallest (3 is found for the most fragile liquids. The 
heterogeneity of time scales suggests possible heterogene- 
ity in space. 

The spatial heterogeneity implied by the mosaic pic- 
ture has received strong support from recent computer 
simulations ^9 20 1 and laboratory experiments, espe- 
cially direct measurements using 4-D NMR and op- 
tical hole burning techniques. In the mosaic picture, 
different regions of the supercooled liquid will relax in dif- 
ferent ways depending on how stable the local structure 
is, but for time scales much longer than the median, the 
system will behave homogeneously since the neighboring 
elements of a mosaic cell will likely also have reconfigured. 
In this picture, despite the presence of some dynamical 
averaging, the response function can be viewed to a good 
approximation as arising from a relaxation time distribu- 
tion P(r) 



^P(r)dT. 



(1) 



This will resemble a KWW formula with the /3 parameter 
determined by the explicit form of P{t), largely by the 
breadth of the distribution. 

According to the random first order transition theory 
|^,^j|jr^ in supercooled liquids the relaxation of an in- 
dividual mosaic element is an activated process. The 
driving force for any local region to escape from one of 
the metastable states predicted by a meanfield free en- 
ergy functional is the configurational entropy of the other 
states to which it might hop. Working against this is a 



1 



cost due to surface energy since the domain wall is not in 
a low free energy configuration. Creating a droplet costs 
free energy that depends on the radius of the droplet. 
Much as in conventional nucleation, one finds 



F{r) = nTscr^ + Ana{r)r'^. 
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(2) 



Here Sc is the configurational entropy density which 
drives random first order transition. A novel feature is 
that the multiplicity of states leads to a renormalization 
of surface tension, cr(r) — o-o(ir)^^^ iQ, where tq is the 
interparticle spacing. The typical free energy barrier is 
determined by the maximum of Eq.(^ as a function of 
r, giving 
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since Sc — Acp(T) , where Ac'p(T) is the specific 
heat jump per unit volume at the transition. The micro- 
scopic theory of ctq gives 
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where Ofi is the square inverse of the Lindemann ratio of 
the maximum vibrational displacement around an aperi- 
odic minimum which is globally stable ulTq = 100. The 
resulting formula for the fragility, D = 32R/Acp fits a 
wide range of substances quite well. (Acp is heat ca- 
pacity jump per mole and R is the gas constant.) The 
typical droplet size in the mosaic corresponding to the 
relaxation barrier turns out to be 
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At Tg , this formula for r 



^ gives a correlation length about 
5 molecular radii (or a few nanometers) |17[] , which is 
consistent with the experimental findings ||21| but much 
larger than the regions envisioned in the Adam-Gibbs ap- 
proach . The configuration of the supercooled liquid 
is separated into domains with average size r* \.Qr^ 
which signals where the free energy Eq. (2) vanishes. 
According to the microscopic theory both and r* (in 
the unit of molecular distance) are nearly universal func- 
tions of the relaxation time at which the response occurs. 
This is because the Lindemann ratio is nearly universal 
for a wide range of substances. Each domain corresponds 
to a minimum of the free energy functional, but these do- 
mains would vary somewhat in size since there are many 
different local minimum states. They will be separated 
by thin mobile sections. The fluctuation of energy of 
each state can also be said to reflect the idea that the 
configurational entropy (at fixed energy!) itself fluctu- 
ates according to the usual Landau formula Q 



(6) 



where = ^i^r^^ is the volume of the average do- 
main. The fluctuation in configuration entropy, the driv- 
ing force, results in a corresponding variation in free en- 
ergy barriers for each mosaic element and therefore gives 
a distribution of relaxation times. 




FIG. 1. An illustration of the "mosaic structure" of super- 
cooled liquids. The mosaic pieces are not necessarily of the 
same size due to the fluctuations in the driving force, configu- 
rational entropy. The system escapes from a local metastable 
configuration by an activated process equivalent to forming a 
liquid-like droplet inside a mosaic element. For droplets with 
size much smaller than the mosaic (as shown with the small 
circle), the droplet shape and its surface energy cost are well 
described by the infinite system result Eq.(2). For transition 
state droplets that would seem larger than the typical mosaic 
elements (as shown with the large circle), the surface energy 
cost will be much smaller than ar^ since the formation of the 
such droplet will break through boundaries of the preexisting 
domain. In this case, the free energy barrier to form such 
a large droplet will be close to AJq, the most probable bar- 
rier determined by the macroscopic configurational entropy 
density. 
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FIG. 2. Correlation between the 13 value for non-Debye re- 
laxation and square root of fragility D. The points are from a 
wide range of experiments collected in [18]. It is important to 
recognize that in this sort of "meta-analysis" of experiments 
that different investigators use somewhat different ways of 
fitting relaxation data and defining the glass transition time 
scale in the laboratory, leading to a spread of values for the 
measured /3 as reported by Bohmer et al. [18]. For exam- 
ple the j3 value for n-propanol (*) is actually measured for 
a glass transition with a time scale of 10~^ seconds [28], 4 
orders of magnitude smaller than the usual glass transition 
relaxation time, 100 seconds, leading to a larger (3. The pure 
Si02 number seems anomalous and deserves careful remea- 
surement since it differs considerably from the Ge02 result. 
The dashed line is obtained assuming a Gaussian free energy 
barrier distribution. It is quite linear. This approximation 
tends to overestimate the width of the distribution resulting 
in a smaller (3. Results from the more realistic cutoff dis- 
tribution taking into account the constraints imposed by the 
mosaic structures yield the solid line. Neither theoretical re- 
sult contains any fitting parameters, owing to the universal 
character of the surface energy costs, obtained by microscopic 
estimates of the domain size of Tg. 

To find the [3 value, we must construct an explicit dis- 
tribution of free energy barriers based on such dynami- 
cally fluctuating mosaic structure. To set the stage, let 
us first see what happens if we assume the distribution 
of relaxation times to be Gaussian resulting from small 
fluctuations. This calculation parallels one carried out 
by Ediger to infer the domain size in supercooled liquids 
psf . This approximation already yields a qualitatively 
(but not quantitatively!) correct correlation of (3 with 
D. The time correlation function can be rewritten as 
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Here t{AF) = tqc^^^''"^. If P{AF) is Gausssian, the 
relaxation function is not precisely a stretched exponen- 
tial but is well fltted by one with a /3 value given by 



/3 = [1 + {5AF/kBTf 



(8) 



where 5AF is the width of the Gaussian distribution 
P{AF) = 1 

(6), we find 
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" — ^ Therefore we have 
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This formula should be valid for a range of substances 
again because of the universal nature of the correlation 
volume at a given time scale predicted by the microscopic 
theory of fragility . The typical barrier height at Tg as 
conventionally defined with a relaxation time of 10^ sec- 
onds corresponds to AFq ~ ilksTg thus this estimate 
gives for [3 at the laboratory glass transition tempera- 
ture ~ The result is also shown on Fig. (2) as a 



dashed line. The data in the graph are the [3 values mea- 
sured at Tg for a wide variety of substances. We see that 
this formula suggests correctly that more fragile liquids 
have smaller /3 values near their glass transition. Also we 
see that as the temperature is lowered, the most prob- 
able barrier height AFq{T) increases rapidly, giving a 
smaller value of (3 ultimately vanishing at Tk if available 
time permitted measurement and equilibration at such 
low temperature. Such a temperature dependence was 
found in the detailed experimental study for o-terphenyl 
by Dixon and Nagel While it seems experiment 

agrees that (3 approaches zero around Tk consistent with 
Eq.(9), some theories do not give such a relation. Free 
volume theory gives a lower limit of /3 = | for example 
i n |27| , a value that has been surpassed in experiments 

m 

Although Eq. (9) gives the right trend oi (3 — D corre- 
lation, it is not quantitatively accurate. A more careful 
analysis of the implications of the mosaic structure is 
needed. The mosaic structure of the random first order 
transition theory implies the existence of large correction 
to the Gaussian result since the fluctuations in AF scale 
in the same way as AFq. They are thus of a similar 
magnitude since is of the same size as r* . We shall 
now show that a reasonably realistic distribution easily 
comes out from a simple model of the dynamic mosaic 
structure of supercooled liquids. First, it is clear that 
even if the distribution of configurational entropy were 
precisely Gaussian, the free energy barrier distribution 
would not be. We should use P{AF)dAF = P{sc)dsc 
to get the precise distribution. Second, domains will not 
all be of equal size, instead, there is a distribution of size 
determined again by fluctuations in configurational en- 
tropy as we have ^ (j")^'^- Most important, when a 
certain domain with small size r' already exists, the free 
energy barrier for overturning that domain will be smaller 
than what Eq.(3) predicts since this equation assumes 
spherical transition state droplets. This activated droplet 
will be modified because of the pre-existing boundaries. 
When the fluctuating size within the spherical droplet is 
comparable or bigger than r' these boundary effects will 
limit the size of the barrier. This effect may be roughly 
described by simply assuming there is a "cutoff' in the 
free energy barrier distribution near the most probable 
one. Similarly if a neighboring region has already flipped 
this will make it easier to reconfigure the domain under 
consideration. Both of these effects suggest the real dis- 
tribution in the mosaic can be better approximated in a 
piece-wise manner, 
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where 

J+°° p{AF)dAF 



a normalization coefficient so 
=1 and P/(AF) is a function deter- 
mined by Pf{AF)dAF = P{sc)dsc. Here AFq is the bar- 
rier assuming no ffuctuation in configurational entropy. 
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This cutoff distribution of free energy barriers still yields 
a (j){t) that fits the KWW formula. The resulting p value 
is, however, different. The results from the cutoff distri- 
bution for (3 at Tg are shown in Fig. (2) as solid line. The 
slope of P versus ^/D is considerably increased and /? now 
saturates to 1 for liquids with D > 150, a value charac- 
teristic of the so-called strong liquids. The improvement 
of agreement over the pure Gaussian is significant, but- 
tressing the case for a dynamic mosaic structure. 

Above Ta, the critical temperature predicted by mode 
coupling theory for dynamic slowing, relaxation is domi- 
nated by collision between molecules. Here dynamics can 
be described by the Debye law with a single exponential. 
We see that between Ta and Tr-, the deviation from the 
Debye law will become more significant upon cooling as 
a consequence of increasing heterogeneity just as in the 
Gaussian analysis. For Gaussian distribution, one finds 
(3 w '^Tk'^ temperature close to Tg. In Fig. (3), we 
also plot the temperature dependence of (3 predicted by 
the more accurate cutoff distribution for o-terphenyl and 
compare it with experiments p^ . In computing the T 
dependence we have assumed that (Tq has its low temper- 
ature value, although near Ta, cr will be smaller, leading 
to faster crossover to the exponential relaxation above 
Ta- 
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FIG. 3. The temperature dependence of /3 for the most 
studied fragile glass former, o-terphenyl is shown. The dashed 
line contains the experimental measurements of Dixon and 
Nagel and their extrapolation to lower temperatures [26] . The 
solid line is the prediction made with random first order tran- 
sition theory. Again no adjustable parameters are present in 
the theory. 

We see that the mosaic structure expected by the ran- 
dom first order transition approach to the glass transi- 
tion can quantitatively explains both the trends in the 
relaxation dynamics over a range of substances and the 
temperature dependence of the deviation from Debye be- 
havior. 
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